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Using the random matrix description of open quantum 
chaotic systems we calculate in closed form the universal au- 
tocorrelation function and the probability distribution of the 
total photodissociation cross section in the regime of quantum 
chaos. 
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Recent advances in laser techniques allow precision 
measurements of the photodissociation of polyatomic 
molecules p| . Realistic quantum calculations of the pho- 
todissociation spectra for small polyatomic molecules, 
such as the radicals HO2 and NO2 Pj or the H^" molec- 
ular ion 0] are also becoming available. The photodis- 
sociation describes the breakup of a bound molecule by 
the absorption of photons, and usually proceeds through 
excited intermediate resonance states. These resonance 
states are directly coupled to the continuum states that 
describe possible breakup channels of the molecule into 
fragments, and their properties are expected to play a 
major role in the dissociation process. Resonances are 
characteristic of all kinds of open systems, i.e. systems 
in which motion along some directions is unbound. 

One of the important observables is the total photodis- 
sociation cross-section. Often it exhibits irregular fluctu- 
ations consisting of partly overlapping peaks (e.g. see 0] 
for the molecule HO2) . Similar fluctuation patterns were 
first observed in resonance neutron scattering [0] and are 
typical for many other systems such as heavy ions Q 
and atoms ||. The statistics of the neutron resonance 
data was explained by random matrix theory (RMT) || , 
where the Hamiltonian is assumed to satisfy the under- 
lying fundamental symmetries of the system but is oth- 
erwise random. Originally justified by the complexity 
of the nuclear system, RMT was later understood to be 
applicable in systems whose underlying classical dynam- 
ics is chaotic ||. Such chaotic systems are expected to 
exhibit fluctuations whose statistical properties are uni- 
versal, i.e. common to systems of different physical na- 
ture and depending only on the symmetry class. RMT 
was successfully applied in the description of both bound 
and open quantum chaotic systems 0. Several univer- 
sal features of chaotic scattering were recently calculated 
in closed form; see Ref. || for details. We therefore ex- 
pect that the photodissociation cross-section of chaotic 



systems will also display universal features. Here we de- 
rive closed expressions for the the autocorrelation func- 
tion and distribution of the total photodissociation cross- 
section in open chaotic systems. 

Following the absorption of a photon, the excited 
molecule can dissociate into several channels. A chan- 
nel describes a fragmentation of the system into sev- 
eral (possibly excited) fragments whose relative motion 
is described by a superposition of incoming and outgo- 
ing spherical waves. We assume that at given energy 
E there are M different possible open channels. A dis- 
sociation solution $>c(E) to the Schrodinger equation is 
denned as a solution that satisfies the following boundary 
conditions: an outgoing wave in exactly one open chan- 
nel c, and incoming waves in all channels. At the given 
energy E there are exactly M independent dissociation 
solutions <I> C (c = 1,2, ... ,M). The total cross section 
<t(E) for the molecule in its ground state \g) (or more 
generally any bound state) to absorb a photon of energy 
E — Eg and to dissociate into any of its M open channels 
is given, in the dipole approximation, by 



a(E) = <r £|($ c (£)|/i|<7)| S 



(1) 



Here (l = fx ■ e is the component of dipole moment /z of 
the system along the polarization e of the absorbed light, 
and cr = (2h 2 e c)- 1 (E - E g ). 

To incorporate RMT description into scattering the- 
ory j7]|| it is convenient to divide the Hilbert space of 
the dissociating system into two parts Q: the internal 
"interaction" region, and the external "channel" region 
where the fragments are far enough from each other that 
their interaction can be neglected. Any solution <£>(£?) at 
energy E can then be represented in terms of its compo- 
nents u and tp, inside the interaction region and channel 
region, respectively. Using standard methods of scatter- 
ing theory (see e.g. ^,0) one can relate the M outgoing 
amplitudes of ip (denoted by the M-component vector B) 
to the inside components u by a linear relation u = CB, 
where 



[E — Hi„ 



rWW f ) 1 W 



(2) 



Here Hi n is the Hamiltonian describing the closed interac- 
tion region when it is decoupled from the channel region 
(e.g. by imposing appropriate boundary conditions), and 
the operator W describes the coupling between the two 
parts of the Hilbert space. 
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Assuming that the classical dynamics of the closed in- 
teraction region is fully chaotic, we can replace the ac- 
tual Hamiltonian Hi n by a random N x N matrix taken 
from the Gaussian orthogonal ensemble (GOE) for sys- 
tems with preserved time-reversal invariance and from 
the Gaussian unitary ensemble (GUE) for systems with 
broken time-reversal symmetry. The coupling W is repre- 
sented by an N x M matrix which we consider to be fixed. 
Because of the invariance of the random matrix Hamil- 
tonian Hi n under orthogonal (unitary) transformations, 
the coupling to the channels is essentially characterized 
by only M invariants, the eigenvalues of W^W. These 
eigenvalues can be expressed in term of the transmission 
coefficients T c (0 < T c < 1), defined through the averaged 
5- matrix by T c = 1 - \ (S CC )\ 2 . The limit T c < 1 corre- 
sponds to an almost closed channel c, whereas T c = 1 
corresponds to the limit of perfect coupling between the 
interaction region and the channel c. 

A central assumption in our model is that direct tran- 
sitions from the ground state to the channels induced 
by the transition operator /x are negligible, and thus the 
decay is possible only via the excited resonance levels. 
Using this fact and Eq. (||), the total photodissociation 
cross-section (|l|) can be rewritten in the following optical 
theorem form: 



a(E)cx(g 



oc Im 
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E — TL 



eff 



9 > , (3) 



where Tt e ff — Hi n — iirWW^ is the effective non- 
Hermitean Hamiltonian known to describe open chaotic 
systems. For closed system (W = 0), Eq. (j|) reduces to 
the strength function of fi, whose statistical properties 
were studied in Refs. [pT|Jl2| . The particular resolvent 
form (|^) has the advantage of being suitable for appli- 
cation of Efetov's supermatrix formalism . The de- 
tailed presentation of such a calculation can be found 
in ||, and here we only present the final result for the 
autocorrelation function of the photodissociation cross- 
section defined by 



S(w = nfl/A) 



(a{E -Q/2)a(E + Q/2)) 
RW 2 
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(4) 



with A being the mean level spacing for the Hamiltonian 
Hi n . S(ui) is found to be a sum of two terms S(lo) = 
Si(uj) + 5*2 (w), which for the GOE case are given by 



dX 



dX 1 



dX 



cos [w(AiA 2 



A)](l-A 2 ) 



[A 2 + A 2 + A 2 -2AiA 2 A-l] 2 



M 



</i,a(A,Ai,Aa)n 



(,9c + A) 



where 



c= L 1 [(5c + A 1 A 2 )2-(A 2 -l)(Al-l)]V2 ' 

(5) 

/i(A,A 1; A 2 ) = 2A 2 A 2 -A 2 -A 2 -A 2 + 1; 



/ 2 (A,A 1 ,A 2 ) = (A 1 A 2 -A) 2 . 

The parameters g c are related to the transmission coeffi- 
cients by g c = 2/T c — 1. We note that each of the contri- 
butions Si^w) represents an interesting object by itself; 
S^w) coincides with the autocorrelation function of the 
Wigner time delays (studied in ||,[l4|-|l6| ) , whereas Si(ui) 
is related to the Fourier transform of the "norm leakage" 
out of the interaction region. The latter quantity was 
introduced recently by Savin and Sokolov as a charac- 
teristic of the process of quantum relaxation in chaotic 
systems [ p7[ . 

For the GUE case one finds S(lu) to be 



S(u) 



dX 



dXi 



Ai 



Ai-A 



n 



g c + A 
g c + Ai 



cos [w(Ai — A)] 
(0) 



In the limit of a closed system (T c = for all c), the 
expressions and (|^) reduce to the strength function 
correlators calculated earlier using RMT and the su- 
persymmetry method ]l2| ]. This limit corresponds phys- 
ically to energies E below the threshold for photodisso- 
ciation (bound-to-bound transitions). Simpler analytic 
forms can be obtained in various limits, such as the limit 
of almost closed systems (i.e. all T c 1), and we defer 
the analysis to a future publication. Here we discuss in 
detail only the limit of a large number of equivalent dis- 
sociation channels (M 3> 1, T c = T) and large density 
of resonances p = 1/A, where all resonances are found 
to have the same width T Such an open system is 
characterized by a single parameter MT = 2irpT = k, 
measuring the degree of resonance overlap [p^,|s||l7|] . In 
this limit of homogeneously broadened resonances one 
finds (cf. 0): 



S 2 (w) 



2 k/2 



(3 u? + k 2 /4 
k/2 

to 2 + K 2 /4 



(7) 

duY 3 , p (u) - , (8) 



where y 2 ,^(^) is Dyson's two- level correlation function 
and f3 = 1,2 for the GOE and GUE cases, respectively. 
For strongly overlapping resonances, i.e. k>1, the func- 
tion S2(u) further simplifies 



1 (k/2) 2 -to 2 
P [w 2 + (k/2) 2 ; 



(9) 



and is only a K _1 -order correction to Si(u>). Generically, 
for many open channels the dominant part of the auto- 
correlation function for lu ~ k is always a Lorentzian. For 
k<1 the tail of S(uj) crosses the u> axis at \ujq\ oc k 1 / 2 3> 
k and then approaches zero from below. We note that in 
numerical simulations of the chaotic photoionization of 
hydrogen atom in external fields the cross-section auto- 
correlation was indeed found to be close to a Lorentzian 
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II . In fact, for M ^> 1 and k 3> 1 it should be possible 
to apply semiclassical considerations and derive the auto- 
correlation function S(cj) from the Gutzwiller trace for- 
mula, as was done for the time-delay correlations fl6U|. 

We emphasize, however, that the autocorrelation func- 
tion for few open channels can differ substantially from 
a Lorentzian. For that purpose, it is instructive to con- 
sider the single-channel case M = 1 where Eq. (||) can 
be reduced to the form: 



S{u)=g h(u) + 



■ 2 

sin uj 



(1 - 2gul 2 (w)) (10) 



with 



. . f°° , COSLUt . . f c 

h(cj)= / dt — — ; J a (w)= / 



sin 
g + l + t 



(11) 



This diverges logarithmically for small ui, dips to a 
minimum below zero ("correlation hole" O]) and then 
exhibits oscillatory decay to zero. As we open the sys- 
tem (i.e. T increases from to 1), the correlation hole 
gradually disappears and the amplitude of oscillations di- 
minishes. Fig. |] shows S(u>) for several values of g. In 
the limit g — » oo, (|Io| ) reduces to S(lj) = 28(u>) — Y2(u) 
p| (shown by the dashed line without the 6- function). 

Parametric correlations jl2) of the photodissociation 
cross-sections can be calculated by incorporating the 
usual factor of exp [-(x 2 /2)/i(A, Ai, A 2 )] (GOE) and 
exp [-a; 2 (A? - A 2 )] (GUE) into the integrand of Eq. (§) 
and Eq. , respectively. 

Another interesting quantity is the distribution func- 
tion of the scaled photodissociation cross-section V(q) = 
(6 (q — a(E)/(a))). To calculate V(q) we make use of 
the observation [ |li^|lO[| that the N x N Hermitean matrix 
CC^ has N — M zero eigenvalues and that its M nonzero 
eigenvalues r c ; c = 1, M coincide with the eigenvalues 
of M x M Wigner-Smith time delay matrix Q = &C. 
Denoting the eigenvectors of CC* that correspond to r c 
by u c , we can rewrite (0) as 



a(E) 



\V C \ 2 ; V c = («?|A|u c ) . 



(12) 



This representation is useful in view of the known statisti- 
cal properties of the time-delay matrix Q f| |l4| , [l5|]2"o"[ . In 
particular, V c , being proportional to the projection of u c 
on the fixed vector jl\g), are independent Gaussian vari- 
ables with the same variance (g\fi 2 \g)/N (V c are real for 
f3 = 1 and complex for (3 = 2). We note that since the 
distribution of the normalized cross-section is indepen- 
dent of the fixed vector (i\g), P{q) coincides with that of 
the quantity p n (E) / (p(E)) , where p n (E) is local density 
of states (LDOS) p n {E) = Im(n\ (E - Heff)' 1 \n) (|n) 
is a fixed vector). Special limiting cases of the LDOS 
distribution were studied in Refs. @Jl|]. 

In the particular case of one open channel the time- 
delay distribution V T (r) is known in the whole crossover 



regime between GOE and GUE symmetries [|15|. In the 
GOE and GUE limits V T {r) simplifies to §: 



P Tj/ j(t) oc t 2+ 2 / d(j)[g + \/g 2 - 1 cos 
Jo 

£ 

2t 



.1/3/2 



(13) 



x exp — — [g + \/ g 2 — 1 cos < 



Using ( |13[) together with the known Gaussian statistics 
of V c and the statistical independence of r c and u c , we 
find from (|l^ ) 



Vp{q) oc 



du 



" u ^, ( q 



q 1 -^ 2 Jo 



[9+ ^fg 2- ! cos <}>f I 2 



q+ [g + \Jg 2 - i cos < 



(14) 
(15) 



For (3 = 1 this expression is equivalent to that obtained 
in 1 19 1 by a different method. 

For the general case of M > 1 open channels, the dis- 
tribution of time-delays r c is known only for the special 
case of perfect transmission (all T c = 1) j20|. In this 
case the LDOS distribution turns out to be intriguingly 
simple Ipf: 



V{q) oc 



(1 + q) W 



(16) 



An alternative way to calculate the LDOS distribution 
is using the supersymmetry approach ]18|] , in which the 
general case of arbitrary transmission coefficients can be 
evaluated for systems with broken time-reversal symme- 
try: 



T{q) = 6(q - 1) + 



47r dq 2 



{2q) 1/2 I dXi 

leff 



1 



(Al - qeff) 1 ' 



i (Ai - A) 

9c + A 
g c + Ai 



(17) 



where q e ff = | (q + q^ 1 )- The integration in ( p7| ) can be 
performed explicitly for an arbitrary set of transmission 
coefficients. For equivalent channels (g c = g) one finds a 
rather simple result: 



V{q) = 



(q 2 + 2gq+l) M + 3 / 2 



[A M (q 2 



1) 



qB 



M 



(18) 



where the coefficients Am , Bm are given by: 

^ = 5[( 9 +i) M+1 -(.«-i)" +, ]S^f 

S M^[( 9+ l)''« + (.,-l)»' + ']2^- 9 4f ■ (19) 
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For the case of perfect transmission g = 1, we repro- 
duce Eq. (^6|). In the opposite limit of an almost closed 
system (g — > oo), we distinguish between two cases: i) the 
regime of isolated resonances (g — > oo at fixed M), and 
(ii) the regime of homogeneously broadened resonances 
(g — > oo and M — * oo in such a way that M/g = At/2 is 
constant). 

In case (i) we distinguish three different regimes for 
q. In the most important intermediate regime X/g <C 
q -C g/M, the distribution is independent of the number 
of open channels M: V{q) oc <7 _1 ' 2 g -3 ' 2 , whereas the 
distribution is A/-dependent for very large and very small 
values of q: V{q) oc g M q~ M+2 for q > g/M, and V(q) oc 
gMqM-i £ Qr y ^ l/q , I n case (ii) one can use the Stirling 
formula in Eq. ( [L8[) to derive the distribution found in 
Refs. [H. 



Except for the one-channel case (14), we are not yet 
able to derive a general formula for the GOE case (i.e. an 
equation analogous to ([l8|)). However, we argue that in 
the limit (i) of sharply non-overlapping resonances (g 3> 
M ~ 1) the correct formula can be obtained by replacing 
in our GUE formula M — ► M/2. The "superuniversal" 
law V{q) oc <? _1 ' 2 g -3 ' 2 is similar to that found earlier for 
the time-delay distribution ||[l5j, and is expected to be 
a generic feature of weakly open chaotic systems. 

In conclusion, we have calculated in closed form the 
universal autocorrelation function and the probability 
distribution of the total photodissociation cross-section 
in the regime of quantum chaos. Our main assumptions 
are the applicability of RMT for describing the closed 
counterpart of the open quantum chaotic system, and 
the absence of direct field-induced transitions from the 
ground state to the continuum. In fact, our results hold 
for the bound-to-continuum strength function of an arbi- 
trary transition operator and for any bound initial state, 
as long as the excited resonance states are fully chaotic. 
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FIG. 1. The autocorrelation function S(ui) of the pho- 
todissociation cross-section vs. lu for g — > 00 (closed system, 
dashed line), g — 20, 5, 2 (solid lines) and g = 1 (perfect trans- 
mission, dot-dashed line). As the system opens (g decreases) 
the minimum gets shallower. The inset shows a magnification 
of the oscillatory behavior of the decaying correlation. 
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